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Temperature dependencies of the impurity magnetic susceptibility, entropy, and 
heat capacity have been obtained by the method of numerical renormalization group 
and exact diagonahzation for the Kondo model with peaks in the electron density of 
states near the Fermi energy (in particular, with logarithmic Van Hove singularities). 
It is shown that these quantities can be negative. A new effect has been predicted 
(which, in principle, can be observed experimentally), namely, the decrease in the 
magnetic susceptibility and heat capacity of a nonmagnetic sample upon the addition 
of magnetic impurities into it. 

PACS numbers: 71.15.Dx, 75.20.Hr 

The Kondo problem is one of most intriguing and important problems in the physics of 
condensed state. It was initially introduced jl| for explaining the minimum of resistance 
in metallic alloys because of scattering of conduction electrons on a magnetic impurity and 
was later extended to other cases. The Kondo effect is key for explaining the behavior of 
heavy- fermion systems [2], the anomalous electronic properties of metallic glasses at low 
temperatures the quantum dots [4|, and many other correlated electronic systems. 

The practically all results concerning this model relate to the case of a band that is flat 
and wide in comparison with the exchange integral J. In particular, this directly concerns the 
analytical solution obtained in terms of the Bethe ansatz However, for real substances 
the electron density of states can have features that cannot be assumed to be wide in 
comparison with J and, therefore, the approximation of a wide fiat zone can be inadequate 
in them. As we show below, the specific features of the electron structure of surroundings can 
manifest themselves in a qualitatively different behavior of the thermodynamic properties of 
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an impurity in such systems. 

Let us write down the Hamiltonian of the Kondo model as follows: 



k(7 



kk' 



As usually j^, let us write down Jkk' = akO^k'J-, where I]k|Q^A;|^ = 1, then extract the 
expression for the local density of states p{uj) = X^k |<^A;P<^(^ ~ ^k), and study the influence 
of its features on the properties of the impurity. 

Following Wilson p we use a unitary transformation to pass from the operators Ck to 
the operators /„: 



Hsd — -J [S^/qi/oi + S /o\/oi + Sz (/o|/oT - /dj/oi 

oo 

"I" 5^ ^^nfnafna + 7" (^fnafn+lcr + fn+laft 



(2) 
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and reduce ([T]) to a semiinfinite chain depicted in FiglH which is called the "Wilson chain" . 
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FIG. 1: Representation of the Kondo model in the form of a semiinfinite Wilson chain. 



ONE-LEVEL BAND 

Let us begin consideration from a very simple case, where the band consists of a single 
level with an energy cq = 0, which corresponds to the case of 70 ^ in Fig{T] This problem 
is easily solved analytically, and we obtain for the partition function of the entire system 
Ztot and for the single level of a free electron Zf^ee- 

Z,,, = 4 + Se^'/' + e-'^'/' , (3) 
where f3 = 1/T, and T is the temperature. 
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The impurity magnetic susceptibility Ximp is defined as a difference in the susceptibihties 
of the entire system, Xtot, and of the system from which the impurity is removed, Xfree- 



Ximp Xtot Xiri 



(4) 



In our case the total magnetic susceptibility of the entire system (we assumed that gfiB = 1 
and /cb = 1) is 



Tx 



1 



tot 



1 + 2e^'^/2 



and the susceptibility of free itinerant electrons is 
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FIG. 2: Magnetic susceptibilities and the impurity entropy for a band consisting of a single level; 
J = -0.3. 

In the case of a negative antiferromagnetic exchange J (at which the Kondo effect does 
take place), the quantity Tximp intersects the abscissa axis with decreasing temperature and 
approaches —1/8 as T (Fig. [2]). This result is understandable: at low temperatures, 
there is formed a singlet with a zero susceptibility in the spinelectronic system; subtracting 
the electronic susceptibility from it, we obtain —1/8. The temperature enters into all 
formulas only in the combination Jf3, therefore, any characteristic temperature will here 
depend on J according to a linear law. 

Apart from the susceptibility (jl]), the so-called local magnetic susceptibility xioc is fre- 
quently introduced as well: 

Xloc(T) = J {S,iT)S^)dT , (7) 
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where {Sz{t)Sz) = ;|Tr e~^^e'^^ S^e''^^ ■ This is the susceptibihty of a single impurity 
in a magnetic field that acts locally, only on this impurity; its magnitude, therefore, can 
hardly be measured experimentally, in contrast to Ximp- We thus obtain 



loc 



^tot 



. 1 BJI2 - 
1 + -e^-^'^ + - 



(8) 
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As we see (Fig. [2]) Ximp and xioc behave quite differently. In principle, this possibility 
was mentioned in Tj, where the reason for the difference is connected exactly with the 
inconstancy of p{uj) and where is asserted that the difference between them disappears for 
a flat band of half-width D in the limit of D — * oo. However, to the best of our knowledge, 
the fact that Ximp can become negative never has been indicated. 

For the impurity entropy, we have 

^in.p = 5tot - ^free = -^{?>{- J / 2)e^'^-^ '^^ + {?>.]/ 2)e~^'^^-^ /'^] + In 1^^^ (9) 

-^tot -^free 

The impurity entropy grows monotonically from — In 4 (at T = 0) to In 2 (at T = oo) (see 
inset in Figl2]). This can easily be understood even without calculations. Indeed, at a zero 
temperature the entropy is equal to the logarithm of the multiplicity of the degeneracy of 
the ground state, while at high temperatures it is equal to the logarithm of the number of all 
states. At T = 0, the entropy of the entire system is Stot = In 1 = (nondegenerate ground 
state) and the entropy of the electronic part is S'^ee = ln4 (four states with an identical 
energy); therefore, S'imp(T = 0) = — ln4. At T = cxo, the entropy of the entire system is 
Slot = ln8 (in the system consisting of an impurity and an electronic site, there are eight 
states); therefore, Sim_p{T = oo) = In 8 — In 4 = In 2. 

As we see, this simple system reveals features that are absent in the well-investigated 
case of a wide flat band, namely, a negative impurity susceptibility and a negative entropy. 
In the general case, we can speak of two limiting cases: wide flat band and inflnitely narrow 
inflnitely high band (single-level band). All the other possible bands are "intermediate" 
between these two, and it is quite probable that some features of the single-level system will 
be preserved also for these "intermediate" cases, which will be considered below in some 
detail. 
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BAND OF A FINITE WIDTH 

Let us examine a symmetrical band of finite width without singularities in an interval 
[-D,D]. 

Analytical Results for J = — oo 

In the case of J = — oo, the Kondo model becomes considerably simpler for calculations. 
When calculating the eigenenergies of the impurity term in (j2]), it proves that the lowest 
energy state is the singlet, and all excited states are separated from it by a gap of about \ J\ 
and, therefore, they can be neglected. Since the matrix elements of the operators /o and /q 
on this singlet are equal to zero, then the neglect of these excited states reduces simply to 
the elimination of members with flfi and ///o from Hamiltonian ([2]) and the replacement 
of the impurity member by its eigenvalue for the singlet state: 

oo 

Hsd{J = -oo) = const + J2 [^nfLfna + In [fl„fn+la + fl+lafna)] , (10) 

crn=l 

SO that the problem is reduced to the problem of noninteracting particles, whose solution is 
much simpler. For example, to calculate Ximp, it is necessary merely subtract the suscepti- 
bility Xfree of the chaiu beginning from the site from the susceptibility of the chain that 
begins from the site 1. 

It is possible to analytically solve the case of a semielliptical local density of states p{uj) = 



2{D^ - uj^y/^/nD^. It is shown in [2| that in this case all the coefficients are 7„ = D/2 
and e„ = 0. The energies of this chain of sites are equal to En = Dcos[7Tn/{N + 1)], 
where n = 1,2, N. We calculate the total density of states V = dn/dE aX n N/2 (i.e., 
at the Fermi level Ep) and obtain T> = {N + l)/7iD. The Pauli susceptibility at the zero 
temperature is equal to x = \1^{Ef). Therefore, 

Ximp( J =---^T = 0) = -—---^ = (11) 

If we now select D = A/tc, so that the height at the Fermi level be p(0) = 1/2, then we 
obtain Ximp = —1/8. 

As far as the impurity entropy and heat capacity are concerned, by using a low- 
temperature expansion for the ideal Fermi gas known from the textbooks, we obtain 

, ^ 27rA^ 27r(A^+l) 27r , , 
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and, according to the well-known formula C = T^, the same for the impurity heat capacity. 
Thus, at low temperatures, we have S'imp = 7imp^ and Cimp = 'jimpT, where 7imp = — 1§; 
i.e., the impurity heat capacity Cimp = Ctot — Cfree is negative. 

An important role in the Kondo problems is played by the Sommerfeld relation, which is 
here more frequently called the Wilson relation 

Ximp/ Xc 47r kj^ Ximp 



R 



(13) 



7imp/7c HQI^b)^ 7imp 

where the subscript "c" indicates the gas of free electrons. In this case R = 1 (this is 
understandable, since the effective Hamiltonian ffTOj) is the Hamiltonian of the ideal gas), 
while Wilsons calculations for small J gave a value R^ 2. 

The same can be made semianalytically for the case of a rectangular local density of 
states: p{u!) = 1/2 for |ti;| < 1. Assuming in formula (12^ (see Appendix) that A ^ 1, we 
obtain the following coefficients for the Wilson chain for a fiat band: 

n + 1 



In 



(2n + l)(2n + 3) 



. 



(14) 



After this, for each length of chain we can easily numerically determine the electronic 
spectrum and, consequently, also the density of states V and the magnetic susceptibility at 
T = 0. 
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FIG. 3: X'm\p{T — 0) for a flat band at J — — oo as a function of the length of the Wilson chain A'^. 



With increasing chain length A^, the magnetic susceptibility Ximp tends to — l/vr^ (oscil- 
lates near this value, approaching it; see Fig. [3]). Analogously, 7imp — ^ —4/3. This differs 
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significantly from tlie above-calculated values for the semielliptical local density of states in 
spite of the fact that at the Fermi level they both have not only equal values, but also equal 
values of the derivatives. 

Given the entire one-electron spectrum, it is possible to numerically determine Ximp, •S'imp, 
and Cimp for any temperature. The results obtained are given in Fig. HI Note separately that 
the impurity heat capacity Cimp is negative, and has a minimum equal to -0.7 at T/D ^ 0.25. 




T/D 



FIG. 4: Tximp, Simp, and Cimp for a flat band of half- width D at J = — oo. 

Emphasize that all above-mentioned data are the result of strict calculations made with- 
out any approximations. 

NRG Calculations for Finite J 

For finite J, the calculations were performed by the standard method of numerical renor- 
malization group (NRG) (see Appendix). Their results are given in Fig. [51 The local 
susceptibility xioc proved to be nonnegative in all the cases, but for ximp there again are 
observed anomalies. 

At \J\ that are considerably less than D (e.g., at J/D = —0.5 in Fig. [S]), we see an 
ordinary picture familiar from the Wilson works. 

At |J| > D (e.g., at J/D = -100 in Fig. [5]), the picture is complicated. At T > D 
(e.g., at T/D > 10 in Fig. [5]), the behavior observed is quite similar to the case of one 
level in the band (cf. Fig. [2]); this is understandable: D is small, and we have p{uj) in the 
form of a narrow rectangle, similar to the 5-function; therefore, the results must be similar: 
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FIG. 5: Impurity magnetic susceptibility at different J. 

Ximp ~ —1/8 at T < |J| and Ximp ~ 1/4 at T ^ oo. At TjD < 10, the behavior of Ximp is 
very similar to Fig. m This is also understandable: since in this case \J\ ^ T, D, the result 
should be close to the case where J = — oo. 

At I J| ~ D, there is a transition region between these two regimes; the low-temperature 
results for it are presented in the table. Note that Ximp and 7imp become zero at different J. 



TABLE I: Low-temperature values of parameters for a rectangular band of half-width D (the errors 
of determining given values are about 1%) 
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Ximp(r = 0) 


Timp 
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-1.08 


0.938 


-1 


-0.005 


-0.5 


0.13 


-0.9 


0.023 


-0.28 


-1.08 


-0.85 


0.034 


-0.22 


-2.03 
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Why have not negative Xi^p and 7imp been revealed in foregoing calculations? This may 
be related to the fact that the procedure of calculating susceptibility used by Wilson was 
suitable only for the case of small J. Wilson used values \ J\ ~ 0.05, and for the low- 
temperature susceptibility he proposed the formula [6|] 



Ximp 



(0.1032 ±0.0005) 



exp(^ \2Jp\ - 1.5824|2 Jp| + 0( JV)) 



(15) 



(where D{Jp) ~ D), which seems to be invalid at large \ J\. From Fig. [5] it is also seen that 
the universality (retention of the form of the curve Tximp{T) at different J) discovered by 
Wilson takes place only for \ J / D\ <^ 1. 

The negative values of Ximp and 7imp appear at J/D ~ —1. This it is explained by a 
change in the behavior of energy levels (Fig. [6l). The greater the magnitude of \J\, the more 
rarefied the energy spectrum and, consequently, the less the total susceptibility Xtot and 
the total entropy Stot- At some, generally speaking, unequal critical values of the exchange 
integral J, they will become less than the independent-of- J Xfree and S'^ee, respectively. 
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FIG. 6: The lowest energies multiplied by A^/^ (A = 2) for the A^th step of the numerical renor- 
malization group, counted off from the energy of the ground state with Q = Sz = {Q is the total 
charge, i.e., the difference between the numbers of electrons and holes; Sz is the z projection of the 
total spin). 
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Piecewise-Constant p{u;) 



Let us now turn to a more realistic model and examine the idealized case of a piecewise- 
constant local density of states with a rise in the center: 



H 



for \u\ < w. 



(16) 



l/2-Hw r ^11^1 
tOT W < \UJ\ < 1. 



The results obtained (Fig. [7]) show that negative heat capacity and susceptibility can 
occur here as well, and even at smaller values of \J\ in comparison with the fiat zone. The 



0.16 
0.14 
0.12 
0.10 
0.08 
0.06 
0.04 
0.02 
0.00 
-0.02 



J = 


-0.25 


J = 


-0.3 


J = 


-0.35 


— J = 


-0.4 




10-^ 



10'^ 



10 



FIG. 7: Tximp and Cimp for a piecewise-constant density of states: the height at the Fermi level 
H = 1.0; the half-width of the peak w = 0.13168. 



general tendency revealed is as follows: the narrower and the higher the central peak, the 
smaller the values of \J\ at which the negative Ximp and Cimp appear. 

And finally, quite a realistic case. The inset in Fig. [8] depicts the local density of states 
p{u!) for one of the sites of the two-dimensional square disordered Anderson lattice taken 
from [8]. The calculation of the magnetic susceptibility of an impurity located at this lattice 
site again shows the presence of a negative Tximp- 
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FIG. 8: Magnetic susceptibilities for an impurity on a disordered t"wo-dimensional lattice; J = 
-0.35. 

LOGARITHMIC SINGULARITY IN THE DENSITY OF STATES 

As is well-known, the density of states p{uj) in two-dimensional systems has a Van Hove 
logarithmic feature. For a two-dimensional square lattice, we have 



where K is the complete elliptical integral of the first kind, which at small uj behaves as 
K[\/l — u"^) ^ In^. It is evident that for this logarithmic feature the characteristic width 
of the central peak is equal to zero; therefore, in the light of what was said in the foregoing 
section it is possible to expect that the negative values of Ximp and S'imp will appear at any 
finite J. 

The Kondo model with a logarithmic feature in the density of states was studied in 
using the perturbation theory in the parquet approximation. In that work, the following 
estimate was given for the Kondo temperature: 



where z/g is on the order of the inverse bandwidth. This qualitatively differs from the widely 
known estimate of Tk for the flat band 




(17) 



(18) 



Tk ~ e 



l/\2Jp\ 



(19) 
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We verified this assertion by NRG calculations, computing via the formula T^XimpiTx) = 
0.0701 that is standard in this method Q]. Then, we performed interpolation using the least- 
squares method with different trial functions (Fig. [9]), which showed that the results much 
better correspond to ( |T8l) than to ( |T9l) . And even more precise is the formula with a root 
correction: 

^ ^ (20) 



with parameters A ^ 10.08 and B ^ 0.206. 
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FIG. 9: Kondo temperature for a two-dimensional square lattice (interpolation by the least-squares 
method): A and B are the variable parameters. 

Using the procedure described in the preceding section, it is possible to obtain strict 
results for the case J = — oo. After making a suitable replacement of variables, it was 
possible to construct simple linear graphs (Fig. [TOl) . from which it is easy to determine the 
low-temperature behavior of Ximp and of Si^p- 
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imp 



T| lnT|" 
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^llnT 



,(v4^ 0.072, a ^0.77) 
J ,(5^ 1.065,5^0.89). 



From fl221) . we obtain for the heat capacity 



a 
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imp 



(21) 
(22) 

(23) 



|lnT|'5+i' 

The NRG calculations carried out for finite J showed (Fig. [TT!) that nothing uncommon is 
observed in the behavior of T^ioc. However, Tximp first intersects the abscissa axis, then 
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FIG. 10: Low-temperature behavior in the case of J = — oo 



begins to grow very slowly, approaching zero from below. Its low-temperature behavior is 
independent of J and coincides with that calculated earlier for the case of J = —oo\ the 
latter is explained by the fact that at sufficiently low temperatures any value of \ J\ manyfold 
exceeds both the temperature and the width of the infinitely thin logarithmic peak of p{uj). 
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FIG. 11: Tximp and Txioc for a square lattice. 



Thus, a new result is obtained: with decreasing temperature, the impurity susceptibility 
Ximp and the impurity entropy Simp/T tend to — oo by the law — 1/T| lnT|"*^'') (Fig. [12]) rather 
than to +00, as in a conventional paramagnet, or to a constant, as in the conventional Kondo 
regime. In this case, the heat capacity is negative and approaches zero according to a 
logarithmic law rather than to a linear law. 
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FIG. 12: Impurity susceptibility and impurity entropy on a square lattice, J = —0.4 

CONCLUSIONS 

In conclusion, let us emphasize the basic result of this work: there was predicted a 
new effect (which, in principle, can be observed in experiment), namely, a decrease in the 
magnetic susceptibility and in the heat capacity of a nonmagnetic sample upon the addition 
to it of magnetic impurities. It must manifest itself in samples in which there are exist 
sharp peaks in the density of one-electron states near the Fermi energy, whose characteristic 
widths do not exceed the magnitude J of the exchange interaction of an electron with an 
impurity. The experimentally observable result must be especially large for the magnetic 
susceptibility in samples with a logarithmic feature in the density of states at the Fermi 
level. 



APPENDIX. DETAILS OF THE METHOD OF NUMERICAL 
RENORMALIZATION GROUP (NRG) 

Here, we briefly describe some moments important for understanding the essence of the 



NRG method (this method has been described in detail in 6j, 10|), as well as some new 
details that were proposed by the author of this work. 

At first, the impurity model with a Hamiltonian of type ([I]) is reduced to a semiinfinite 
chain with a Hamiltonian of type ([2]). The logarithmic discretization and constructing a 
Wilson chain are conducted numerically (except for the calculations of the flat band, which 
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can be made analytically). 

The renormalization-group procedure starts from the solution of the isolated-impurity 
problem (sites "imp" and eo in Fig. [T]). At the first step, there is added a first conducting 
electronic site ei, and there is constructed and diagonalized a Hamiltonian matrix on this 
Hilbert space (with a fourfold higher dimensionality). This procedure is repeated multiply. 
However, since the dimensionality of Hilbert space grows as 4^ (A^ is the order number of 
an iteration), it is not possible to store all eigenstates during the computation. Therefore, 
it is necessary to preserve after each iteration only states with the lowest energies. For 
the model in question, if we limit ourselves by a certain maximum number of stored states 
(determined by the possibilities of the computer at hand), it is necessary, beginning from a 
certain iteration, to leave on the order of 1/4 of states at each step. 

Unfortunately, the direct application of this diagram meets with failure, since the distur- 
bance introduced by the rejection of the high-lying states proves to be too large. Wilson 
found a method of overcoming this difficulty; it reduces to the artificial introduction of an 
exponential damping of matrix elements 7„, which decreases the coupling between the pre- 
served and rejected states, thus decreasing the infiuence of the rejected states. To this end, 
Wilson used a logarithmic discretization of the conduction band, i.e., replacement in ([1]) 
of an energy range G {DA~~"~^,DA~^], n = 0,1,2,... by a single level with an energy 
equal to the average energy of this interval {D is the half-width of the conduction band, 
A > 1). As a result, after passage from (P) to ([2]) the jumps will have the required damping: 
7n oc A~"/^. For a fiat band, Wilson obtained analytically 

D(l + A-i)(l-A-"-i) 



7r. 



A-"/2 , e„ = . (24) 



2Vl - A-2"- VI - A-2"-3 
The algorithm of obtaining the result for the initial nondiscretized band is as follows: to 
conduct the calculation of the observed impurity value at a given A, changing the number 
of states preserved at iteration, and to make sure that the result is independent of it, i.e., 
it is reliable for a given A; to make this for different A (for example, for A = 3, 2.5, 2, 
1.5) and, using those values for which it was possible to obtain reliable results, to make an 
extrapolation A ^ 1. 

Now, we concern the new details of the method that were suggested by the author. Earlier, 
when numerically passing from ([I]) to ([2]), the researchers attempted to numerically conduct 
a unitary transformation from the operators Ck to the operators /„. For this, it is necessary 
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io solve an infinite set of linear equations, which is impossible; therefore, for example, in 
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m 



the number of equations was artificially limited. We instead used a method (described 



of constructing a Wilson chain via three-diagonalization according to Lanczos 12| . 
which consists in iteratively constructing such a space (the so-called Krylov space) in which 
the operator 

H^ = J2 ^kcLcko- (25) 

k 

has a three-diagonal matrix. Starting from the vector |0) = /oo.|vac) (where /g^ = Skt^fcCjLcr) 
and using the formula 

H^\n) = -fn\n + 1) + en\n) + -f*Jn - 1) (26) 

we consecutively obtain the coefficients 7„ and e„ for ([2]). 

As is indicated in 6|, because of the retention of only part of the energy spectrum at 
the iVth step of the renormalization group procedure, the thermodynamic averages should 
be calculated at a temperature that depends on A: T/v = A~^/^To, where the starting 
temperature Tq is selected more or less arbitrarily. We made an attempt to select the 
last value in a regular way. In a chain of finite dimension, the thermodynamic quantities 
experience evenodd oscillations depending on (see Fig. [3]). It was discovered, that the 
magnitude of these oscillations depends on the choice of Tq and has a clearly pronounced 
minimum; the temperature corresponding to this minimum was selected as starting. 



To finally suppress these evenodd oscillations, we used smoothing according to Euler [13| , 
which reduces to the following. If there is a certain oscillating sequence An, we introduce a 
new sequence A!^\ whose members are equal to arithmetic mean of the adjacent members of 
the initial sequence: Aj;^^ = {A^ + An+i)/2. If necessary, the procedure is repeated: A'^^ = 
{A^'^ + Anli)/2. In particular, this was made in the calculation of Ximp- By designating 
Xn = Ximp(TAr), we obtain 

X^^KV'^nTn+i) = ^Xn + ^Xn+1 
X^^KTn) = \xN-i + \XN + \XN+1 (27) 
X^^HV^A^^Af+l) = \xn-i + \xn + \xn+i + \xn+2 

With the aid of this method, we determined the low- temperature values of Ximp and 7imp, 
which are given in the table [11 
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